#include "Matrix.h"
#include <iostream>

int main(int argc, char *argv[])
{
    Matrix<double> A(2,3);
    A(0,0) = 0;
    A(0,1) = 0.1;
    A(0,2) = 0.2;
    A(1,0) = 1.0;
    A(1,1) = 1.1;
    A(1,2) = 1.2;
    std::cout << "A = \n";
    for(int i =0; i <A.get_n_row(); i++)
    {
        for(int j=0; j< A.get_n_col(); j++)
            std::cout << A(i,j) <<" \t";
        std::cout<<"\n";
    }
    std::valarray<double> p(3);
    p[0] = 1;
    p[1] = 2;
    p[2] = 3;
    p = A * p;
    for(int i =0; i < p.size(); i++)
        std::cout<< p[i] << "\t";
    std::cout<<std::endl;
    Matrix<double> A_T = A.transpose();
    std::cout<<"A^T = \n";
    for(int i =0; i <A_T.get_n_row(); i++)
    {
        for(int j=0; j< A_T.get_n_col(); j++)
            std::cout << A_T(i,j) <<" \t";
        std::cout<<"\n";
    }
    return 0;
}